From damage percolation to crack nucleation through finite size criticality 
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We present a unified theory of fracture in disordered brittle media that reconciles apparently 
conflicting results reported in the literature. Our renormalization group based approach yields a 
phase diagram in which the percolation fixed point, expected for infinite disorder, is unstable for 
finite disorder and flows to a zero-disorder nucleation- type fixed point, thus showing that fracture 
has mixed first order and continuous character. In a region of intermediate disorder and finite system 
sizes, we predict a crossover with mean-field avalanche scaling. We discuss intriguing connections 
to other phenomena where critical scaling is only observed in finite size systems and disappears in 
the thermodynamic limit. 
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Brittle fracture in disordered media intertwines two 
phenomena that seldom coexist, namely, nucleation and 
critical fluctuations. The usual dichotomy of thought be- 
tween nucleated and continuous transitions makes the 
study of fracture interesting. Even more intriguing is 
the fact that crack nucleation happens at zero stress in 
the thermodynamic limit: smaller is stronger and larger 
is weaker. This makes the existence of critical fluctua- 
tion in the form of clusters and avalanches of all sizes 
even more mysterious. What kind of critical point gov- 
erns a phase transition that happens at zero applied field 
(stress) in the thermodynamic limit, and what is the uni- 
versality class of such a transition? How do self-similar 
clusters, extremely rough crack surfaces, and scale in- 
variant avalanches ultimately give rise to sharp cracks 
and localized growth? These questions have been ad- 
dressed previously via a host of different theories, such as 
those based on percolation and multifractals [lH^] , spin- 
odal modes and mean- field criticality jEj], and classical 
nucleation In this Letter, we present a theoreti- 

cal framework based on the renormalization group and 
crossover scaling that unifies the seemingly disparate de- 
scriptions of fracture into one consistent framework. 

Fracture in disordered media is the result of a complex 
interplay between quenched heterogeneities and long- 
range stress fields leading to diffuse damage through- 
out the sample, and local stress concentration favoring 
the formation of sharp localized cracks. The self-afflne 
morphology of crack s 11 (J, the power-law statistics of 
avalanche precursors [Ill4l4j, and the scale dependence 
of the failure strength distribution 15- 1?}, all result from 
this competition. Disordered fracture can be understood 
in the limit of infinitesimal as well as infinite disorder. 
Infinitesimal disorder means perfect crystalline material 
with just a few isolated defects (say a missing atom or 
a micro-crack). In this limit, fracture statistics can be 
understood as a nucleation type first order phase transi- 
tion 6-91. In the limit of infinite disorder, stress concen- 



tration becomes irrelevant and fractureprogresses via un- 
corrected percolation-like damage [l|, l2| • This mapping 
to percolation theory becomes rigorously valid when the 
disorder distribution is not normalizable (or very broad, 
in the language of multifractals) [lj. The situation is 
more interesting at intermediate disorder, where unlike 
typical first order transitions, crack nucleation is pre- 
ceded by avalanches with power-law distributions and 
mean- field exponents [H|, [l8l - [2o| . sometimes interpreted 
as a signature of a spinodal point (BJ. Our renormaliza- 
tion group based theory unifies the above descriptions 
into a single phase diagram. 

We use a 2D fuse network to model disordered brit- 
tle materials. A description of the disordered fuse net- 
work model that we study can be found in any number 
of references 1211, El, EH 13 . Briefly, we consider a peri- 
odic network of fuses arranged in a square lattice of size 
L tilted by 45° (the so-called 'diamond lattice', figure 
la). Each fuse is assigned a quenched current threshold 
from a common distribution with a cumulative distribu- 
tion function F(-). If the current through a fuse exceeds 
its threshold, then the fuse is burned and is removed 
from the network i.e., its conductance is set to zero. The 
current through the network is ramped quasi-statically, 
and fuses are burned one at a time until the network 
becomes non-conducting, at which point the network is 
said to be fractured. We assign thresholds between and 
1, specifically we take F(x) = x& , f3 > 0. This form of 
distribution of thresholds serves as model for a generic 
distribution with a power-law tail at the origin, and has 
been studied widely [||-[5| . In this model the limit f3 — >> 
corresponds to infinite disorder, while the limit f3 — » oo 
corresponds to infinitesimal disorder. Figure [1] shows a 
schematic of an undamaged fuse network (QJO, and re- 
alizations of fractured networks for various values of the 
parameter /?. Notice how the damage looks percolation- 
like for small f3 while a single crack appears for large /?. 

We begin by arguing that crack-tip stress concentra- 
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(a)A fuse network (b)/3 = 0.03 
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(c)/3 = 0.5 (d)/3 = 3.0 

FIG. 1: Fuse network model, a). Schematic of a fuse net- 
work. Periodic boundary conditions are used in the horizontal 
direction, b-d). Fractured sample for various values of the 
parameter /?; the spanning cluster (or crack) is colored red. 
There is a smooth crossover from percolation-like behavior for 
small f3 to nucleated cracks at large /3. 



tion is a relevant perturbation to the infinite disorder 
percolation critical point. Our assertion implies that 
percolation-like behavior is a finite-size crossover effect. 
This is consistent with recent results showing that even 
an arbitrarily small cutoff in the threshold distribution 
(at f3 = 0) leads to a crossover away from percolation at 
large system sizes [23|. We calculate self-consistent up- 
per and lower bounds for the stress and damage fraction 
at failure, and show that all of these quantities vanish 
in the limit of large L. This will establish that percola- 
tion cannot be the dominant behavior for large L, since 
percolation demands that the damage fraction be finite. 
Let (7 f , (j)f be the stress and the damage fraction at fail- 
ure, respectively. The lower bound on both quantities is 
trivially equal to 0. The upper bound is obtained self- 
consistently. Let us assume that <pf < 0^ <C 1, where 
0yT is an upper bound on 0/, and similarly af < <C 1. 
Let, if possible, the damage be percolation-like, so that 
0j = F{a~f) 30] . The stress at the tip of a crack of 

length / (lattice units) is given by au p (l) ~ a~f (1 + aV7), 
where a is a lattice dependent constant. Thus, the length 
of a critical crack at a given stress and damage fraction is 
lcr{o~f) °° l/( cr /) 2( ^ 2 + h.o.t. The probability that a crit- 
ical crack forms at a given lattice site is at least F(a^) lcr . 
Since there are L 2 sites in the lattice, the probability of 



1 such crack appearing on the entire lattice is at least 
L 2 F(a~f) lcr |6[. At the failure stress this probability is 1, 

thus <jj can be obtained by solving L 2 F(a^) lcr< ' (J f ^ = 1. 
It can be proved that the solution cr^(L) — >• as L — >• oo, 
thus, 0+ = F(cr+ ) H> 0. 

We have established that percolation is unstable to nu- 
cleation, however, the crossover length is expected to be 
rather large. The reason for this effect is that 4>~f(L) 
decays very slowly with L. This slow decay is expected 
to also manifest itself in the form of large finite size ef- 
fects. The rate of decay obviously depends on F(-), for 
F(x) = x? one can show that </>+(L) - (/3/2 logL)^/ 2 . 
More sophisticated estimates that account for stress con- 
centration during the growth of the critical crack, as op- 
posed to percolation-like growth assumed here, yield sim- 
ilar results. The convergence becomes extremely slow as 
f3 approaches 0, meaning that percolation threshold will 
be reached before nucleation of the critical crack for in- 
creasing larger system sizes. This is consistent with the 
previous studies that found that the fuse network can 
be mapped onto a percolation problem in the limit of 
f3 — >• [2|. However, one should note the subtle point that 
order of limits matters since percolation is ultimately un- 
stable to nucleation at any f3. 

The avalanche behavior associated with fracture can 
be understood via a simple argument. The argument 
is valid in the vicinity of the critical point and breaks 
down for very large L. Consider an avalanche that 
starts with a bond breaking at a stress a (<C 1) and 
damage fraction <j> (<C 1). Linear elasticity predicts 
that the change in the stress field due to the break- 
ing of the bond, c(r, a), decays as c(r, a) ~ cr/r 2 (ig- 
noring the dipolar directional dependence), where r is 
distance from the broken bond. The probability that 
a bond at distance r breaks in response to this change 
in stress is approximately given by F'(a)c(r, a). Thus, 
the expected number of bonds that break in response 
to stress change due to one bond breaking is given by 
A ~ rdrF f (a)c(r^a) ~ F'{cr)cr logL. Substituting the 
form F(x) = x 13 gives A(a, (3,L) ~ ficjP log L = f3(j)logL. 
This shows that lirn^o 0, L) =0 (for fixed L), thus 
there are no avalanches for small /3, and the damage is 
percolation-like. For suitable /3 the avalanche progresses 
as a branching process, where breaking of one bond trig- 
gers a few more and so on (A is also known as the branch- 
ing ratio). It is well known that integrated avalanche size 
distribution for such processes is a power-law with expo- 
nent r a = 3/2 + 1 = 5/2; for suitably large L we expect 
the avalanche size distribution to be a power-law with 
exponents consistent with the mean field value of 5/2 [H[ . 
Finally, for very large L (or /?), the system flows away 
from the critical point and the avalanches get cutoff due 
to nucleation effects. 

All the ideas discussed so far can be encapsulated 
neatly in the form of crossover scaling functions. The 
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scaling form for the cluster size distribution can be de- 
rived by using ideas of scale invariance. Let G(zi, . . . , z n ) 
be a scale invariant function, then by definition, G(-) 
should remain invariant under a rescaling by a factor 6, 
i.e. G(zi, . . . , Zn) = b a °G(zib ai , . . . , z n b an ) for some con- 
stants OLi. Taking b = 1 + e and solving up to first order 
in e gives the general form of a scale invariant function 
as G(z 1 ,...,z n ) = z^ ao/ai g(z 2 z^ a2/a \. . . ,z n Zi an/ai ), 
where the universal scaling function, £/(■), and the crit- 
ical exponents, c^/ai, are characteristic of the critical 
point [3l|, [1^. The variables Z{ represent directions in 
parameter space near the critical point. The directions 
with OLi > belong to the relevant parameters and those 
with oti < to irrelevant parameters. We treat /3,1/Lto 
be a relevant parameters, and let u be the leading irrel- 
evant parameter (the largest of the negative c^). Thus, 
ignoring all irrelevant variables but the leading one, the 
scale invariant distribution of cluster sizes can be written 
as 

P c (s\P,L) =s~ Tc T c (pi}'"* isL-V^f , u L~ A f^ 

where the subscript c denotes variables associated with 
the clusters. We use the subscript / (for fracture) to 
distinguish the critical exponents from their counterparts 
in percolation theory. We know that in the limit of j3 —> 
(at fixed L) the cluster size distribution should reduce to 
distribution of percolation clusters at the critical point, 
thus we can deduce three critical exponent combination, 
namely r c = 187/91 = 2.0549, a c v f = 48/91 = 0.5275 
and Af/i/f = 72/48 = 1.5 [HI . Even thought the clusters 
created in fracture are loopless [23| , the static properties 
of loopless percolation are identical to usual percolation, 
thus the use of percolation critical exponents is valid [26| . 
The moments of the cluster size distribution should scale 
as (taking a Taylor expansion in uL~ A f/ Uf for large L) 

{S U C ) = iCn+l-rc)/^!/, (^jc^ L l/u f ^ L -A f /u f 

where JZ (-), /Q(-), n = 2,3..., are universal scaling 
functions [32|. From a data fitting perspective, it is easier 
to deal with the moments (as opposed to the distribution 
function) because J7^(-), A^(-) are functions of just one 
scaling variable. The functions for the avalanche size 
distribution are completely analogous, 

P a {s\p,L) = s~ r ^ a ^L^f.sL-^f.uL-^/^ , 

where r a is expected to be close to its mean field value 
of 5/2. 

We have done numerical simulations to verify our the- 
oretical predictions. We did extensive statistical sam- 
pling of systems of size up to L = 128 and /3 between 
0.03 and 8. In order to fit the data to the scaling predic- 
tions we use the following functional forms for the scaling 
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(a) Avalanche size distribution (b) Cluster size distribution 




(c) Scaling collapse of (s^) (d) Scaling collapse of (sj?) 



FIG. 2: Scaling theory of fracture, a). The avalanche 
size distribution shows a power-law consistent with the mean 
field exponent of 5/2 for moderate f3 (= 0.5 at L = 128). As 
expected, the power-law is distorted for much smaller or larger 
f5. b). The cluster size distribution shows a power- law that is 
consistent with the exponent predicted by percolation theory 
(= 187/91). The power- law cutoff becomes smaller as one 
moves away from the critical point, c, d). The scaling forms 
fit the data well, confirming the predictions of the scaling 
theory. Higher moments of the distributions fit the scaling 
forms as well (not shown here). Notice the significant finite- 
size effects as the data gets closer to the L = oo curve with 
increasing system sizes. 



functions for the moments of the cluster size distribution 
(with y c (x) = (log a; - /i c )/a c ) 

i—m 

J*(x) = a , n evf(y c (x)) + e -^(*)) 2 £ (y c (x)) , 

i—m 

JC c n (x) = ai, n erf (y c (x)) + e -(*(*» a £ B ? n Hi (y c (x)) , 



where /x c , a c , a , n , «i,n, A^ n , B c i n are fitting param- 
eters, erf(-) is the error function, and Hi(-) is the i th 
Hermite polynomial. We use the first three Hermite 
polynomials in the expansion, i.e., m = 3. The cor- 
responding forms for the avalanches are (with y a (x) = 
(log a; - Ha)/oL a ) 



j:(x)=e-^^ 2 J2 A lnHi(y a (x)), 

i=m 

K a n {x) = e -to«(-)) 2 Bt, n H t (y a (x)) . 



i=0 



The forms of the scaling functions are chosen so that 
they have the correct asymptotic behavior. As discussed 
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FIG. 3: Phase diagram for brittle fracture in disor- 
dered media [4]. Disorder decreases along the /3 axis; 
nucleation governs the behavior for small disorder or long 
length scales. Percolation is characteristic of the large dis- 
order regime, while the crossover region exhibits interesting 
critical behavior in the form of scale free distributions of 
avalanche sizes. The topology of the fractured samples evolves 
from percolation-like damage for large disorder to well de- 
fined sharp cracks in the nucleated regime. The phase bound- 
aries are quantitatively somewhat arbitrary, and are set at the 
value of the scaling variable fiL 1 ^* at which the second mo- 
ment of the avalanche size become half of its peak value (for 
the avalanche phase); the boundary of the percolation phase 
is found analogously. 



previously, we know that lim 
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since there are no avalanches for very small j3 (at 
fixed L) and at very large L (at fixed /?). On the 
other hand we know lim^i/^Q J^L 1 ^) = C, for 
some constant C (according to percolation theory) and 
lim^ L i/ ! y / ^ oo J^(fiL x l Vf ) = since there are no clusters 
in the nucleation dominated regime away from the critical 
point. The forms used here satisfy all these requirements. 

Figure [2] shows the size distributions as well as fits to 
the scaling forms. It is evident that the data is consistent 
with the scaling theory. Based on joint fits for n = 2, 3 
(n = 3 not shown in figure [2]) we estimate the follow- 
ing values of the critical exponents: Vf = 1.56 ± 0.30, 
a a = 0.47 ± 0.15, A f = 2.35 ± 1.50, a c = 0.34 ± 0.08. 
For the fits shown in figure [2] the exponent r a is held 
at its mean- field value of 5/2, while unbiased fits yield 
r a = 2.45 ± 0.25. The scaling exponent r c , and the ex- 
ponent combinations <r c ^/, Af/vf are held at their the- 
oretical values of 187/91 and 48/91, 72/48, respectively. 
The statistical error bars are much smaller than the error 
bars reported here. We have estimated the error bars due 
to systematic errors by using a variety of techniques such 
as varying the number of terms in the scaling functions, 
trying different fitting forms, varying the critical range 
for the fits, varying the error bars on the data over a 
reasonable range, etc. Figure [3] shows the phase diagram 



that emerges from our analysis [33] . In the (5 — l/L space, 
curves along which the scaling variable fiL x l Vf attains a 
critical value demarcate the boundary between qualita- 
tively different behavior. Note that the exact position of 
the boundaries is somewhat arbitrary, since this is a not 
an abrupt (first order) transition; however, the diagram 
is qualitatively accurate. 

The critical phenomena associated with fracture has 
several intriguing characteristics. Firstly, the scaling 
function associated with the avalanches has a singular- 
ity at 0, linw/^Q J£(/3L 1 / l/ f) = 0, that subdues the 
avalanche behavior as the critical point is approached. 
Secondly, there is no point in the phase diagram (ex- 
cept for the j3 = limit) that shows any critical phe- 
nomena in the limit of L — >• oo. Thus, scale invar iance 
itself becomes a finite-size effect; perhaps this phenom- 
ena should be named finite-sized criticality. Finally, it 
is rather remarkable that the critical phenomena (typi- 
cally associated with continuous phase transitions) gives 
way to nucleation (a first order transition) in the limit 
of long length scales! Thus, fracture has mixed first or- 
der and continuous transition character. Transitions of 
mixed first order and continuous character have become 
somewhat of a theme in the past decade or so. Re- 
cently we noted that the Mott transition and dielectric 
breakdown have a mixed character [27j; similar findings 
have been reported in a variety of fields such as jamming 
transitions, rigidity percolation [28], and phase-separated 
manganites [29| . 

In conclusion, we have presented a scaling theory of 
fracture that builds on renormalization group ideas and 
unifies several disparate results in the field. Our the- 
ory shows that percolation-like behavior as well as the 
scale invariant precursor avalanches leading to fracture 
are finite-size effects. We show that on long length scales 
brittle fracture is always nucleated. We hope that our 
analysis will pave the way for a deeper understanding of 
the many mysteries associated with the phenomenon of 
fracture. 
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